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Nonrelativistic bound states are studied using an effective field theory. Large logarithms in the effective theory 
can be summed using the velocity renormalization group. For QED, one can determine the structure of the 
leading and next-to-leading order series for the energy, and compute corrections up to order a 8 In 3 a, which are 
relevant for the present comparison between theory and experiment. For QCD, one can compute the velocity 
renormalization group improved quark potentials. Using these to compute the renormalization group improved it 
production cross-section near threshold gives a result with scale uncertainties of 2%, a factor of 10 smaller than 
existing fixed order calculations. 



1. INTRODUCTION 

Nonrelativistic bound states in QED and QCD 
provide an interesting and highly nontrivial prob- 
lem to which effective field theory methods can be 
applied jl]|2|. The QCD bound states we will con- 
sider are heavy QQ states such as it bound states 
or the T system. In QED, the classic examples 
are Hydrogen, muonium (/i + e~), and positron- 
ium. Each of these systems has three important 
scales, m the fermion mass, mv the fermion mo- 
mentum, and mv 2 , the fermion energy. (For Hy- 
drogen and muonium, m is the electron mass or 
the reduced mass of the two particles.) The ve- 
locity v is of order the coupling constant (ev s or 
a), and we will only consider the case v -C 1, 
mi) 2 3> Aqcd so that nonperturbative effects are 
small. 

Multiscale problems with widely separated 
scales are well suited for study using effective field 
theories. For example, if the problem has the 
scales mi ^> m 2 ^> m 3 . . ., one first starts with 
the theory above mi, and matches to an effec- 
tive theory below m± in which only modes with 
masses much smaller than mi are retained. The 
effective theory is then scaled using the renormal- 
ization group to the next scale m 2 . At this point, 
particles with mass m 2 are integrated out to con- 
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struct a new effective theory, and so on. The 
complicated multiscale computations of the orig- 
inal theory are reduced to a number of simpler 
single scale computations of matching and run- 
ning in the effective theory. The effective theory 
method also allows one to sum logarithms of the 
ratio of mass scales lnmj/mi+i using the renor- 
malization group evolution between m^ and m I+ i. 

The goal is to correctly separate the scale m, 
mv and mv 2 for nonrelativistic bound state prob- 
lems using an effective field theory, and to sum 
large logarithms using the renormalization group. 
The large logarithms in this case arc lnp/m, 
hiE/m and hip/ ' E which are proportional to hiv, 
and lead to In a contributions to bound state en- 
ergies. Furthermore, for QCD, the effective the- 
ory also determines the scale of the strong cou- 
pling constant, i.e. whether one should use a s (m), 
a s (mv) or a s (mv 2 ). The nonrelativistic effective 
theory, NRQCD/NRQED, has been studied ex- 
tensively in the past @-|lJ]. What is new is the 
precise formulation of the effective theory, and the 
way in which the renormalization group is scaling 
is implemented. 

The results presented here will be applied to 
the study of it production in the threshold region. 
There is a large ratio of scales, m t ~ 175 GeV, 
m t v ~ 26 GeV and m t v 2 ~ 4 GeV, where 
v ~ 0.15 is the typical velocity in the nonrela- 
tivistic bound state. Clearly a s lav is not small, 
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and summing logarithms is important in this case. 

The results arc also useful in QED. While 
a In a is small, it is important to compute 
to high orders because the experiments have 
high precision. The Hydrogen Lamb shift of 
1057.845 MHz is known to an accuracy of 
9 KHz [Q, the Hydrogen hyperfinc splitting is 
measured to be 1420.405 751 766 7(9) MHz fl9), 
and the muonium hyperfinc splitting is 
4463.302 776(55) MHz @. The binding en- 
ergy of Hydrogen, m e a 2 /2 is 2 x 10 10 MHz, so 
the experimental error in the Lamb shift of 10 
ppm is a part in 10 12 of the binding energy. 
We will be able to compute corrections of order 
m e a 8 In 3 a/ (An) 2 ~ 5 KHz to the Lamb shift, 
which are relevant for the present comparison 
between theory and experiment. [The counting 
of 47T factors for bound states is a little different 
than the conventional counting |21| . Potential 
loops give powers of a whereas soft and ultrasoft 
loops give powers of a/(4vr).] 

A detailed comparison of theory and experi- 
ment for QED can be found in Refs. &Q. 



2. NEW RESULTS 

There are many interesting new results that 
have been obtained for QED and QCD |l5|j2| 
|2l| . For QED, one finds a universal description 
of In a terms. A single renormalization group 
equation gives the Lamb shift, hyperfinc splitting 
and decay widths for Hydrogen, muonium and 
positronium. The renormalization group method 
allows us to compute for the first time the a 8 In 3 a 
Lamb shift in positronium and the a 8 In 3 a Lamb 
shift in Hydrogen and muonium including recoil 
corrections. It also resolves a controversy in the 
literature about the a 8 In 3 a Hydrogen Lamb shift 
in the limit m p — > oo. 

The renormalization group method allows one 
to understand the structure of the QED pertur- 
bation series, and why the In a corrections termi- 
nate. The leading order series has a single term 
that contributes at order a 5 In a to the energy, 
and the next-to-leading order series terminates 
after three terms, a 6 In a, a 7 In 2 a, and a 8 In 3 a. 
One also finds some infinite series of terms in 
QED, but they have the form (a 3 In 2 a)™, rather 



Figure 1. A potential gauge boson exchange. The 
typical momentum and energy transfered are mv 
and mv 2 . 



than (a In a)™. 

In QCD, one is able to sort out the scales for a s , 
and decide whether the strong coupling is a s (m), 
a s (mv), or a s (mv 2 ). We also obtain the renor- 
malization group improved computations of the 
bound states potentials in QCD. There arc nu- 
merous applications of these results, and I will 
show an example of the dramatic improvement 
one obtains for the tt production cross-section 
near threshold p7j , p8| ]. 

3. THE PROBLEM 

The basic problem can be seen by drawing a 
few Feynman diagrams. A typical gauge boson 
exchange in the t channel such as Fig. |l| has mo- 
mentum transfer of order p ~ mv. A wavef unc- 
tion graph or radiated gauge boson graph such 
as Figs. @ have gauge boson momenta of order 
E ~ m\r . More interesting diagrams such as 
those in Figs. [?] involve gauge bosons with mo- 
menta of order p and order E. In a graph such as 
Fig. ^, the vacuum polarization insertions make 
the effective coupling of the two gluons a s (mv) 
and a s (mv 2 ) respectively. One result which 
should be clear from Fig. || is that graphs can in- 
volve a s {mv) and a s (mv ) simultaneously. We 
will return to this important point later on. 

4. MOMENTUM REGIONS AND DE- 
GREES OF FREEDOM 

The Feynman integrals in the full theory can 
be evaluated using the threshold expansion p9]| . 
The important momentum regions (in Feynman 
gauge) are referred to in the literature as hard 
(E ~ m, p ~ to), potential {E ~ mv 2 , p ~ mv), 
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Figure 2. Graphs containing ultrasoft photons, 
with energy and momentum of order mv 2 . 
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Figure 3. Graphs containing gauge bosons carry- 
ing momentum of order mv and rav 2 . 
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Figure 4. An example of a graph involving both 
a(mv) and a(mv 2 ). 



ultrasoft (E ~ m» 2 , p ~ rau 2 ) and soft (15 ~ raw, 
p ~ mv). The threshold expansion momentum 
regions are often used to describe bound state 
computations; however it is important to note 
that the threshold expansion is not an effective 
field theory. To construct an effective field theory, 
one needs to include only modes that can be on- 
shell. The effective theory therefore has nonrela- 
tivistic fermions (which are potential modes), and 
soft and ultrasoft gauge boson modes. The hard 
fcrmion and gauge boson momentum regions, the 
soft fermion momentum region, and the poten- 
tial gauge boson momentum region do not require 
modes in the effective theory. 

The desired effective theory is valid for ener- 
gies and momenta much smaller than the fermion 
mass m. One can try expanding in powers of 
E/m and p/m as in heavy quark effective the- 
ory, so that the expansion parameter is 1/m. For 
example, the dispersion relation E = y / p 2 +m? 
gives terms in the Lagrangian of the form 

^ (e~^- + ^ + ..) V. (1) 



L 



8m 3 



The lowest order propagator is 1/(E + ie), which 
gives 9(t) in position space. This is the static 
propagator of HQET: fermions propagate forward 
in time, but do not move in space. This propa- 
gator is acceptable for some calculations involv- 
ing heavy quarks. For example, one can com- 
pute the static potential between fixed sources 
using this propagator. However, for it produc- 
tion, the quarks are produced at the same point, 
and they remain at the same point for all time if 
the static propagator is used. This is too singu- 
lar, and the HQET expansion breaks down. In 
general, it is essential for treating nonrelativistic 
bound states that the heavy fermions move. For 
this to occur, the lowest order propagator should 
be 1 / (E — p 2 /2m + ie) , so that E and p 2 /2m are 
of the same order in the effective theory power 
counting. This implies that the 1/m expansion 
cannot be used; instead one must use an expan- 
sion in powers of v, where E and p 2 /2m are both 
of order v 2 |],|. 

The effective theory expansion parameter is the 
velocity v, and formally, a must also be treated 
as order v. Thus order a 2 radiative corrections to 
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the leading term are just as important as order 
v 2 relativistic corrections. The effective theory 
below the scale m has: 

• Nonrelativistic fermions with propagator 

1 

E - p 2 /2m + it 

• Ultrasoft gauge bosons coupled via interac- 
tions that are multipole expanded ||. 

• Potentials V(p,p') for the scattering of an 
incoming Q and Q with momenta p and — p 
to outgoing Q and Q with momenta p' and 

P'- 

• Soft gauge bosons. The importance of in- 
troducing soft fields in the effective theory 
was first pointed out by Griesshammer ]30| . 

The effective theory has two different gauge bo- 
son fields, soft bosons and ultrasoft bosons. This 
does not lead to any double counting if graphs are 
evaluated in dimensional regularization. 

The static theory is not the m — ► oo limit or the 
v — > limit of the effective theory. For this rea- 
son, the static potential and the effective theory 
potential are not equal. 

5. POWER COUNTING 

The power counting parameter of the effective 
theory is the velocity v. If one expands the dis- 
persion relation as in Eq. (]]]), then E and p 2 /2m 
are both of order v 2 , and p 4 /8?n 3 is of order v 4 , 
i.e. of order v 2 relative to the leading term. 

The potential V(p,p') also has an expansion 
in powers of v. The leading term is the Coulomb 
potential, V(p, p') oc a/ |k| 2 , where k = p' — p 
is the momentum transfer. Since momentum is 
of order mv, the Coulomb potential is naively 
of order a/v 2 . However, the potential is a four- 
fermion operator, whereas the kinetic energy is 
a two-fermion operator. This leads to an addi- 
tional factor of v from the power counting factors 
for the fields, so that the Coulomb potential is of 
order a/v in the effective theory. One can then 
determine the power counting for all the other 




Figure 5. An iteration of two potentials in the 
effective theory. 



potentials by comparing with the Coulomb po- 
tential. The hyperfine interaction oc aSi • S2/1T1 2 
is generated by one-photon exchange, and is of 
order v 2 relative to the Coulomb interaction, so 
it is of order av in the power counting, as are the 
spin-orbit, tensor and contact (Darwin) interac- 
tions. At one-loop, there are also potentials that 
are proportional to odd-powers of k. The first 
such potential is proportional to a 2 / |k|, and is of 
order a 2 v in the power counting. 

A loop graph such as Fig. [5] of the time-ordered 
product of two potentials of order a ai v bl and 
a a-2 V b 2 i s f or d er a ai+a 2v b 1 +b 2 _ Qnc can now see 

that the static potential differs from the m — > 00 
or v — ► limit of the effective theory potentials. 
For example, the loop graph of Fig. || with one 
l/(m|k|) and one Coulomb potential is of order 
a 2 v° x a/v = a 3 /v, and is of the same order 
in v as the Coulomb potential. The two particle 
intermediate state propagator 1/(E — p 2 /2m) = 
2m/ (2mE — p 2 ) produces a factor of m in the nu- 
merator, that cancels the 1/m at the vertex. In 
the static theory, the l/(m |k|) potential is set to 
zero before the loop integration, so that the graph 
of Fig. U is not present in the static theory. As a 
result, the NRQCD potential Q differs from the 
static potential. 



6. MATCHING CONDITIONS 

The method of calculating matching conditions 
is the same as in any effective theory. One com- 
putes the graphs in the full theory at the scale 
/x = to, and subtracts the corresponding graphs 
in the effective theory. The graph in Fig. |^ gives 
the matching condition for the fermion potential. 



5 




Figure 6. Tree-level matching for the potential. 




Figure 7. One-loop matching for the potential 
in QED. The first line gives examples of the full 
theory graphs. The second line gives examples of 
effective theory graphs: an iteration of two po- 
tentials, and a soft photon graph. The difference 
of the two sets of graphs gives the one-loop cor- 
rection to the potential. 

The full theory amplitude, 

(p-p') 2 

is expanded in powers of p,p', to give the po- 
tential in the effective theory. At one-loop, the 
difference of the full theory and effective theory 
graphs in Fig. |^ give the one-loop corrections to 
the matching potential. The only difference at 
this stage between Hydrogen and positronium is 
that there arc annihilation contributions to the 
positronium potential from graphs such as Fig. |[ 
The graphs can have an imaginary part, that give 
the positronium decay width. 




Figure 8. Annihilation contributions to the 
positronium potentials. The second graph has an 
imaginary part. [There is also a one- loop crossed 
box in the annihilation channel.] 




Figure 9. An ultraviolet divergent two- loop graph 
involving the iteration of three potentials. 



7. RENORMALIZATION GROUP EVO- 
LUTION 

The nonrelativistic bound state system has 
three important mass scales, m, mv and mv 2 . 

7.1. Two-stage running 

The conventional method of implementing the 
renormalization group is as follows 

• Start at /i = m 

• Scale \i from m to mv 

• Integrate out the soft modes at mv 

• Scale fj, from mv to mv 2 

This is referred to as the two-stage method, 
because there are two-stages of renormalization 
group evolution. Consider a loop graph involv- 
ing time-ordered products of potentials, such as 
Fig. ^. This graph contains a logarithm of the 
form In \J mE / ji. When /j is set to mv, this log- 
arithm has the from In \J E/mv 2 , and is small. 
Thus the logarithms in the graph are summed by 
renormalization group evolution of \x from m to 
mv. 
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Figure 10. One- loop ultrasoft photon renormal- 
ization of the potential. 



The graph in Fig. |lOj involving an ultrasoft pho- 
ton exchange contains a logarithm of the form 
hxE/fi. The fi in this ultrasoft graph is scaled all 
the way down (in two stages) to tow 2 , at which 
point the logarithm is In E/mv 2 , and also small. 

However, this two-stage method of implement- 
ing the renormalization group turns out to be in- 
correct for nonrelativistic bound states. The rea- 
son is that the scales mv and mv 2 are correlated — 
one cannot be varied independently of the other. 
Instead one needs to use an alternative one-stage 
scaling procedure. 

7.2. One-stage running 

In one stage running, one introduces two differ- 
ent n parameters, /15 and fijj JnJ. In dimensional 
regularization in 4 — 2e dimensions, the soft pho- 
ton coupling is multiplied by fi s , the ultrasoft 
photon coupling by \x\j, and the potentials by 
/iff. Note that this is only possible because we 
have two different photon fields to represent the 
soft and ultrasoft photons in the effective theory. 
Then 

• Set us = mu, fijj = mv 2 

• Start at v = 1 and scale to v = v. 

This procedure is referred to as the velocity renor- 
malization group, because one runs in veloc- 
ity v rather than momentum |ll[. The loga- 
rithms in Figs. ^| and [l^ are now In \fmE/mv 
and In E/mv 2 , which arc minimized when v = v. 
Thus this method also minimizes logarithms in 
the diagrams, and sums them by renormalization 
group evolution. 

The difference between the two renormalization 
group methods can be seen in Fig. HJ. In 
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Figure 11. Paths in the (fiu,fj,s) plane for one- 
stage and two-stage running. 



two-stage running, there is only a single (i, so 
that us = fJ-u = M, and they are lowered together 
from to to mv. At this point, the soft modes are 
integrated out, and for the ultrasoft modes is 
lowered to mv 2 . The integration path in Fig. |ll| 
is along the lower edges of the triangle. In one- 
stage running, the integration path is along the 
diagonal. It is convenient to define two anoma- 
lous dimensions, 7s and "fu by taking the deriva- 
tives of Green's functions with respect to In us 
and In respectively. One can show by explicit 
calculation that 

• The two paths give different answers. The 
integration is path dependent because V x 
7^0. 

• One-stage running using the velocity renor- 
malization group agrees with explicit QED 
calculations at order a 3 In 2 a, a 7 In 2 a and 
a 8 In 3 a. 

The moral is that for nonrelativistic bound states, 
one should run in velocity rather than momen- 
tum. 

The difference between the two integration 
methods can be made more precise. In the two- 
stage method, one first integrates 7s + from 
li = to to /i = mv, and then integrates "fu from 
/x = mv to fi = mv 2 . In the one-stage method, 
one integrates 7s + 2ju (since In fj,jj runs twice as 
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fast as In fig) from v = 1 to v = v. If the anoma- 
lous dimensions arc constant, the two methods 
give 



(75 + 7l/)ln^+ 7 t/ln^ 



two- stage 



(75 + 27(7 ) hi v one-stage 

and agree with each other. However, in general 
anomalous dimensions are not constant, but can 
depend on coupling constants Vi, that themselves 
run. As a result, one finds that the lnv terms 
agree, but the higher order terms differ. For ex- 
ample, consider a In 2 v term that depends on the 
product of 75 and 7c/. For two-stage running, 
the contribution is proportional to ^slu + ®7u = 
Islu from the two pieces of the path. For one- 
stage running, the contribution is 7s (27^), which 
differs by a factor of two. Similarly, a ^fslu hi 3 v 
contribution differs by a factor of four, and so on. 



Table 1 

Numerical values for the it singlet potentials. 
The values at v = 1 are the matching values at 
/x = m t . The values at v = v are the veloc- 
ity renormalization group improved values, where 
v = 0.14 has been used. 



Coefficient 


v = 1 


V = V 




-1.81 


-2.47 


mU ( k s) 


-0.36 


-0.03 


m 2 U [ r s) 


-1.81 


-1.49 


m 2 U { 2 s) 





0.63 


m 2 u[ s) 


0.60 


0.53 


m 2 U [ l ] 


0.15 


0.16 


m 2 U ( t s) 


2.71 


3.11 



8. RUNNING POTENTIALS 

The running potential V(p, p') has an expan- 
sion 

v(p, p') = + + + v& + ... (3) 

where is of order v n in the velocity power 
counting. The first three terms in the expansion 
have the form 

k 2 ' 



1/(0) = £k 

V (1) =U 2 + U a S 2 



(4) 



U r (p 2 + P' 2 ) 



2k 2 



i\J A ■ (p' x p) 
k 2 

-UJtrx ■ cr 2 



3 k • CTi k • er 2 
k 2 

where - 1/m, and ~ 1/m 2 . In QCD, 
each of the coefficients can be written as U — > 
U^l ® 1 + U^T A <g> T A , where 1 and T A /f A 
are color matrices acting on the quark/antiquark 
lines. The anomalous dimensions for the coeffi- 
cients U c -Ut have been computed, and the de- 
tails are given in Refs. H-& The renormal- 
ization group improved static potential was com- 
puted in Ref. pjj. An important point to note 



is that graphs can involve both soft and ultrasoft 
gluons, so that the anomalous dimensions involve 
both a s (mv) and a s {mv 2 ). As an example, the 



running of U. 



(i) 



is given by 
14Ci 



/mv\ 



3 

3/3 



a s (mis)a s (m) In ( 



a s (m) In 



a s (mv) 
a s (mv 2 ) 



(5) 



where C x = 2/9 for QCD. Note that Eq. (§) de- 
pends on a s (m). a s (mv) 1 and a s {mv 2 ). The run- 
ning coefficients in the singlet channel (C/( s ) = 
— CfU^) for it production are presented in 
Table [l]. The renormalization group improved co- 
efficients U r and U2 differ significantly from their 
matching values, because they depend on the ul- 
trasoft scale through a s (mv 2 ). The other coeffi- 
cients only have a soft anomalous dimension, and 
do not run as much. 

The renormalization group improved poten- 
tials can be used to calculate the renormaliza- 
tion group improved cross-section for it produc- 
tion in the threshold region. Fig. |l2| shows a 
sample fixed order calculation of i?, the ratio of 
a{e+e- -> tt)/a{e + e~ -> f|. The scale 

uncertainty is of order 20%. The renormalization 
group improved version of the results is shown 
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as 




0.0 t — i 1 1 1 1 1 1 1 1 1 1 1 . 1 1 — J 

346 347 34S 349 350 351 352 353 354 

Figure 12. Fixed order computation of ft produc- 
tion near threshold. The curves are LO (dotted), 
NLO (dashed) and NNLO (solid.) Uses the 15 
mass-scheme [[33] [35| 




o.o t — i 1 1 1 1 1 1 1 1 1 1 1 — . 1 1 — J 

346 347 34S 349 350 351 352 353 354 

Figure 13. Renormalization group improved com- 
putation of it production near threshold. The 
curves are LL (dotted), NLL (dashed) and NNLL 
(solid.). Uses the 15 mass-scheme [^3|~^H| 



in Fig. 13. There is a dramatic reduction in the 
scale uncertainty, which is now around 2%, as well 
as an improvement in convergence for the normal- 
ization. The small theoretical uncertainty means 
that an accurate measurement of the cross-section 
can be used to study new physics. For example, 
a standard model Higgs boson of mass around 
115 GeV changes the cross-section by ~ 5%, and 
is measurable. 



9. QED 

The velocity renormalization group method 
gives very interesting and important results when 
applied to QED jl5| . The basic potentials we will 
need for QED are summarized in Table |2|. The 
last column gives the contribution to the bound 
state energy levels due to the given potential. The 
fourth column gives the order of a given poten- 
tial, treating v as order a. Since the Coulomb 
potential is of order unity, one finds the obvi- 
ous result that the Coulomb potential must be 
summed to all orders, and cannot be treated as 
a perturbation. The potentials yW, V^, 

are first generated at tree-level, and are of order 
a, whereas the potentials V^°\ V^- 2 \ are 
first generated at one- loop, and are of order a 2 . 

The bound state energy levels can be deter- 



mined to order a 4 by computing the matrix ele- 
ments of and between Coulomb wave- 
functions. Time-ordered products of two po- 
tentials, such as T[y(°V (0) ], T[y(°V (1) ] and 



Table 2 

Table of potentials for QED. The first column is 
the potential, the second gives typical terms in 
the potential, the third gives the power counting 
in a and v, the fourth gives the order in the v 
counting scheme when v ~ a, and the fifth gives 
the contribution of the potential to the bound 
state energy. 







Power Counting 


Order 


E 


y(-l) 


a 
k 2 


a 

V 


1 


a 2 


y(o) 


a 
m k 


a 2 


a 2 


a 4 


y« 


a aS 2 
m 2 1 rn 2 


av 


a 2 


a 4 


y(2) 


a k 

m 3 


2 2 

a v 


a 4 


a 6 


y(3) 


ak 2 

m 4 


av 3 


a 4 


a 6 
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T [yWyW] first contribute at order a 6 . In prin- 
ciple, to obtain the energy levels to order a 4 , one 
also needs the one- and two-loop matching cor- 
rections to the Coulomb potential. However, such 
corrections vanish in QED. As a result, the first 
correction to the order a 2 binding energy is of 
order a 4 , and is given by the matrix element of 
y(o) There are no order a 3 corrections to 

the energy levels in QED. 

Define the leading and next-to-leading order 
anomalous dimensions of a potential to be the 
anomalous dimension from graphs at one and two 
higher orders in a than the potential itself. For 
V^ -1 ) and V^ 1 ' which are of order a, the leading 
order anomalous dimension is of order a 2 , and the 
next-to-leading order anomalous dimension is of 
order a 3 . For V^ ' which is of order a 2 , the lead- 
ing order anomalous dimension is of order a 3 , and 
the next-to-leading order anomalous dimension is 
of order a 4 . Since different terms in the poten- 
tial arc of different orders in a, the terms leading 
and next-to-leading order are not related to the 
number of loops. 

Integrating the renormalization group equa- 
tions for and V™ using the leading order 
anomalous dimension gives a series of the form 



[ 1 + a In a ■ 



a 2 In 2 a 



; ln d 



which contributes 



a 4 (1 + a In a + a 2 In 2 a + a 3 In 3 a + . . .) 

to the energy. Integrating the next-to-leading or- 
der anomalous dimensions gives 

a 4 a (l + a In a + a 2 In 2 a + a 3 In 3 a + . . .) 

terms in the energy. The next-to-next-to-leading 
anomalous dimension gives 

a 4 a 2 (l + a In a + a 2 In 2 a + a 3 In 3 a + . . .) , 

terms in the energy, which are the same order as 
those obtained by using the leading order anoma- 
lous dimension for the and potentials 
which first contribute at order a 6 . Thus one can 
compute the 

a 5 In a a 6 In 2 a a 7 In 3 a 
a 6 In a a 7 In 2 a a 8 In 3 a 

series in the energy using 7lo, 7nlo for V*- 0,1 -*. 




Figure 14. One-loop running of in QED. 
The first graph has a soft photon, and the other 
graphs have ultrasoft photons and a potential. 



10. LEADING ORDER 

The Coulomb potential and do not run in 
QED at leading and next-to-leading order, so one 
is left with the running of V^'. The anomalous 
dimensions are evaluated for a particle of mass mi 
and charge — e interacting with a second particle 
of mass TO2 and charge Ze. Evaluating the graphs 
in Fig. [l4] gives 



dU 2 
dv 



UZ 2 a 2 



2a 



3mi7772 37T 



1 

7771 



z 

7772 



(6) 



where the first term is the soft contribution from 
Fig. ^4] a and the second is the ultrasoft contribu- 
tion from Fig. |l4|b,c. Note that the ultrasoft con- 
tribution has been multiplied by two, since the 
anomalous dimension for the velocity renormal- 
ization group is 75 + 2777. The other coefficients 
in (U r , U s , Ua, Ut) have zero anomalous di- 
mension at this order. 

Since the Coulomb potential and a do not run 
in QED, one can combine the two terms, 



dU 2 
v—r- = 7oU c 
dv 

which defines 
2a ( 1 



^° 3TT V m l ' 477717772 ' 777^ 



(7) 



(8) 



7o is a constant in QED since a does not run. 
Integrating Eq. (^) gives 



U 2 (v) = U 2 (l)+~( Q U c lnv, 



(9) 



where U2 is evaluated at v = v = a. Since 70 
is a constant, Uiiy) only has a ln^ term, and 
terms of the form In" is, with 77 > 1 vanish. As a 
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result, the leading order energy series Eq. (|9|) ter- 
minates after a single term, so one has an a 5 In a 
contribution to the energy, but the a 6 In 2 a, etc. 
terms vanish. At low orders, the absence of terms 
other than a 5 In a in the leading order series has 
been noticed before by an explicit examination of 
Fcynman graphs. This is the first general proof 
that all the terms beyond a 5 In a in the leading 
order series vanish 

The matrix element of U 2 gives the energy shift 



AE = (U 2 {v)) 

= 7o t/ c lnH^(0)| 2 
8Z 4 a 5 m|j 
37rn 3 

where we have used 
|#)f=^ 



1 



4miTO2 m. 



(10) 

— * I In Za, 



(11) 



for the nS state, and ttir is the reduced mass. 
This is the famous a 5 In a correction to the Lamb 
shift first computed by Bethe, including all recoil 
corrections. 

11. NEXT-TO-LEADING ORDER 

At next-to-lcading order, the anomalous di- 
mension for V^ 1 ' is 



dU- 



2+s 



dv 



Pccc Ul + Pcc2 U 2 C (U 2+s + U r ) 



NLO 



+p c22 U c (u 2 +s + 2U 2+8 U r + hj? - 5C/ 2 S 2 

+Pck U c U k + p k 2 U k (U 2+s + U r /2) 

+Pcs Uc (u 3 + U 3s S 2 + lu rk ] , (12) 



where U 2 + s = U 2 + U S S 2 , and the coefficients are 



Pccc 

Pc21 
Peel 
Pc3 



1 

647r 2 ymf 

,2 



4tt 2 ' 



1 



1 



8w 2 \rrii m f 2 
2m R 



Pck 
Pk2 



2tt 2 
2mu 



rrii 



ml 



(13) 



The anomalous dimension Eq. (|12j) can be inte- 
grated by substituting the leading order running, 
Eq. (^) for the coefficients on the right-hand side. 
Since only U 2 runs at leading order, the right 
hand side has at most a In 2 v, so that the in- 
tegral has at most a In 3 v term. This implies that 
the next-to-leading order series Eq. (^|) terminates 
after the first three terms, a 6 In a, a 7 In 2 a, and 
a 8 In 3 a. 

11.1. In 3 a 

The only term that contributes to the In 3 a cor- 
rection is the t/f term of Eq. ([l2"|). Integrating 
gives a contribution to U 2 (y) of the form 

^Pc 22 f/ C 3 (l)ln 3 ^ (14) 

which is spin- independent, and has no imagi- 
nary part. There is no contribution to the decay 
width or hyperfine splitting at this order. The 
Lamb shift at this order is obtained by multiply- 
ing Eq. (|l4| ) by the matrix clement of the unit 
operator, |t/;(0)| 2 , to give 

R ■ In 3 (Za) 



AE = 



27n 2 n 3 
1 



Z 2 



Am\m 2 



(15) 



which is approximately 8 KHz for the 2P-2S 
Lamb shift in Hydrogen. Substituting Z = 1 and 
mi = m 2 = m e gives the a 8 In 3 a Lamb shift for 
positronium 



AE 



3ra e Q 8 In 3 a 



(16) 



87r 2 n 3 ' 

The positronium Lamb shift is a new result, 
as are the recoil terms in the Hydrogen Lamb 
shift. In the limit m\/m 2 — > 0, the Hydrogen 
Lamb shift has been computed previously by sev- 
eral groups. There is an analytic computation 
by Karshenboim f36j | and a numerical compu- 
tation by Goidenko et al. |5?J that agree with 
our result. There are also numerical computa- 
tions by Malampalli and Sapirstein [B8|, and by 
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Figure 15. Four-loop diagram that contributes to 



the a In a Lamb shift 



Yerokhin which agree with each other, but 
disagree with the other results. Recently, there 
has been a computation by Pachucki |4(j| that 
agrees with our result. Yerokhin [^4j has em- 
phasized that the complete a 8 In 3 a Lamb shift 
might not be contained in the loop-aftcr-loop cal- 
culations of Rcfs. 

The other calculations rely on extracting 
the logarithm from four-loop diagrams such as 
Fig. |l5j. The velocity renormalization group fac- 
tors the graph into the product of a two-loop 
anomalous dimension p C 22, and the square of a 
one-loop anomalous dimension 7q. 



11.2. 



In 2 a 



The In a contribution to the hyperfine split- 
ting and decay widths is given by the U2{v) con- 
tribution 



7op c22 C/ c 2 (l) [C/ 2 (1) + C/ S (1)S 2 



In 2 v 



(17) 



The spin-dependent term gives the a 1 In 2 a hy- 
perfine splitting 



HFS = - 



9mim27m 3 
1 Z 



Z 2 



4m 17712 



(18) 



where pi are the magnetic moments normalized 
to unity for a Dirac fermion. Our result agrees 
with previous calculations |36|]42[ | . Substituting 
the matching values for U s (l) for positronium 
(which differs from Hydrogen because of annihi- 
lation contributions), one finds the positronium 
hyperfine splitting 



Ps HFS = - 



7m e 

87T77, 3 



7 1 2 

\ m a, 



(19) 



which agrees with a recent computation of Mcl- 
nikov and Yclkhovsky 43 1. The imaginary 
parts of the matching coefficients give the decay 
widths 1 36 1, 



Ar 



7o Pc22 U c (l) 2 In 2 v 



3 a 3 ln 2 a, (20) 



for both ortho- and para-positronium. 
11.3. In a- 

The In a contributions to the decay width arise 
from 



U 2+s 



pc22U c (U 2 +s+2U r ) 



+Pcc2U 2 + P2kU k 

which give 



In v 



(21) 



f " fa***- 



7S 2 

6 



— 2 a 2 lna, 



(22) 



so that 

1 / ortho 

/Ar N 



a 



In n 



= -2a 2 In a. 



(23) 



These agree with existing results |44 4 



12. CONCLUSIONS 

The methods presented here give a systematic 
way of separating scales in nonrelativistic bound 
state problems. All large logarithms are summed 
using the velocity renormalization group. The 
method provides a universal description of QED 
logarithms. The agreement with known results 
at order a 5 In a, a 6 In a, a 7 In 2 a, and a 8 In 3 a is 
a highly non-trivial check of the formalism. In 
QED, one finds that the leading order series ter- 
minates after one term, and the next-to-lcading 
order series terminates after three terms. In ad- 
dition, the method resolves a controversy about 
the a 8 In 3 a Lamb shift for Hydrogen, and gives 
the first calculation of the a 8 In 3 a energy shift 
for positronium. 
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In QCD, one can distinguish a s (mv) and 
a s {mv 2 ), and both can appear simultaneously in 
the same anomalous dimension. The renormal- 
ization group improved potentials can be used to 
compute it production, and reduce the scale un- 
certainties by a factor of ten. 

The velocity renormalization group should also 
be applicable to other problems with correlated 
scales. In the bound state problem, one can gen- 
erate the scale mv in loop graphs from the scale m 
and mv 2 , mv 



j 2 . Similar effects can 



occur at finite temperature, where one has the 
scales T, gT and g 2 T, and some of the ideas de- 
scribed here might be applicable to that problem 
as well. 

This work was supported in part by the De- 
partment of Energy under grant DOE-FG03- 
97ER40546 and by NSERC of Canada. 
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